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Abstract 

The Lanczos process has been analyticahy and exactly 
carried out for the spin 1/2 isotropic XY chain in the 
thermodynamic Hmit, yielding a form for the Lanczos 
coefficient 0^{s). This coefficient has a monotonic vari- 
ation for real positive s and confirms a general theorem 
on the ground state properties of extensive Many-body 
Systems. The Taylor expansion of the coefficient about 
s = has a finite radius of convergence, and ground state 
estimates based on a finite truncation of this are shown 
to be asymptotic. 



PACS: 05.30.-d, ll.15.Tk, 71.10.Pm, 75.10.Jm 
The Lanczos Algorithm is one of the few reliable 
general methods for computing ground state and some 
excited state properties of interacting quantum Many- 
body Systems. It has been traditionally employed as a 
numerical technique on small finite systems, with atten- 
dant round-off error problems, although the main ob- 
stacle to its further development is the rapid growth of 
the number of basis states with system size. We wish 
to demonstrate that the Lanczos Method can be used 
exactly and analytically for quantum Many-body Sys- 
tems, and in addition where the thermodynamic limit 
can also be treated exactly. We demonstrate this for a 
member of an exactly solvable class of such systems, those 
falling into the universality class of the two-dimensional 
Ising Model, namely the isotropic spin 1/2 XY Model. 
In fact we have provided another notion of solvability, 
both in the mechanics of constructing the solution and 
in terms of its general applicability. In addition we state 
some general theorems arising from the Analytic Lanczos 
Method Formalism, which apply to both integrable and 
non-integrable models, where in the latter case a trunca- 
tion in an expansion necessarily takes place. Given the 
exact results for our model the effects of this truncation 
can be examined in minute detail, and we present such a 
convergence study. 



In our convention the Lanczos process generates an 
orthonormal basis \ijjn) from a suitable trial state ji/jo) 
with the recurrence |jl|, 0| 

iJlV'n) = I3n\^n-l) + a.n\lpn) + /3„+l|V'ri+l) • (1) 

All of our attention focuses on the Lanczos coefficients so 
defined, a„ and /3^, rather than the basis states and we 
perform these steps exactly. It was found recently^, |j 
that the coefficients arising from an extensive Many-body 
System satisfy a confluence property an{N)/N —f a{s) 
and /3^{N)/N'^ as n, ^ oo where the scaled 

Lanczos iteration number is s = n/N. With such a de- 
scription, and employing theorems on the extremal zeros 
of Orthogonal Polynomials, a number of exact general 
theorems for the ground state energy |3|, ground state 
averages and gapsQ, |^ of an arbitrary extensive Many- 
body System were found. For example the ground state 
energy density is given by 



eo = inf [a(s) 

s>0 



(2) 



For the N site chain, our model Hamiltonian is taken 



to be 



N 
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(3) 
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where the label i index the sites, and where we only con- 
sider the c-cyclic problem, as the thermodynamic limit 
will be taken. This model, being at the isotropic point, 
is at a critical point of an anisotropic family of models, 
and consequently has vanishing excited state gaps. The 
trial state was taken to be the classical Neel state even 
though this state is in fact, as far from the true ground 
state as one could possibly be but its simplicity allows 
for ease of moment computation, certainly compared to 
alternatives. So as a "natural" choice we proceed by only 
considering this choice. In the spin picture the Neel state 
lAf) will be represented by | t|tJ, . . . titi)) where the 
quantisation axis is in the z-direction. If we denote the 
vacuum state for the quasiparticles by |0) then the Neel 



1 



state can be expressed in terms of the quasiparticle op- 
erators 77J by 



1-^) = n ^^10) • 

qez ^ ^ 



(4) 



Here the spectrum of ahowed momentum values have 
been divided up into the set Z = [—n, —n/2) U (0,7r/2) 
and its complement Z, and there is a bijection from Z to 
Z defined by g = —q—ir if q < and q = —q+n if g > 0. 
Basically the Neel state is created from the quasiparti- 
cles vacuum state by pairs of excitations whose momenta 
sum to ±7r. 

One first constructs the cumulant generating func- 
tion from its Taylor series expansion about t = 0, 



(AA|e*^|AA)=cxp|5]^.„-Ue 



(5) 



and the exact generating function F{t) was found to beg 

F{t) = — dq\ncosh{tAn) , (6) 

Jo 

where the quasiparticle energy is A, = | cos(g)|. So from 
the Taylor series expansion about t = we have the 
cumulant s 



22«-l(2n)! 



An {n\Y' 



n > 1 



(7) 



where B2n are the standard Bernoulli numbers. 

Our moment generating function is also essentially 
the characteristic function of the underlying measure 
defining the Lanczos Process, so a Fourier inversion in- 
tegral yields the weight function w{e) in which the only 
parameter is the system size N. In what follows we will 
only need the leading order contribution as iV ^ cx) so 
we perform a stationary phase approximation yielding 



f 1 r/^ 

lnu;(e) = TV < ^e-l- - y (i(?lncosh(^A, 



(8) 



where ^(e) is the stationary point determined from 

dqAq tanh(^Ag) , 



(9) 



and e = E/N the energy density. 

A key link to complete our analysis is a very recent 
result from the theory of random Matrix Ensembles j§] in 
which the Lanczos coefficients are determined from the 
weight function to leading order in n, TV ^ 00 by the 
supplementary condition 







a+20 



de 



u'ie) 



2p ^4/32 - ' 



(10) 



and the normalisation condition 
Q+2/3 
= — / de- 



27r 



a-2/3 ^4/3' - 



(11) 



where w(e) = — Inu'(e). This pair of equations arises from 
the continuum Coulomb gas approach to the eigenvalue 
distribution of random Matrix Ensembles, in the solution 
for the density of eigenvalues cr(e), where the u(e) plays 
the role of a confining one-body potential. The explicit 
solution for this charge density can be constructed, but 
we do not require this here. An equivalent derivation [pO[ 
of the above results, which does not involve Matrix En- 
sembles, employs a continuum Coulomb gas approach 
to the distribution of stationary points arising from the 
leading order stationary phase approximation to the Sel- 
berg integral expression for the Hankel determinants. 

In our model a{s) — and our fundamental result 
then is the following pair of implicit equations for 0^{s) 
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2/3 



de- 



■k/2 



(12) 



dq cos(q) tanh(^ cos(q)) 



In Figure |T| we have plotted the dependence of /? on 
s in a form which illustrates the ground state theorem, 
Eq.(||), (eg = — I/tt). Clearly — 2/3(s) approaches eq as 
s — > 00, i.e. there is no minima but an infimum, and 
the monotonicity of this Lanczos coefficient can be easily 
shown from the defining equations 



ds 
d^ 



?0 



(13) 



with 2/3 = e(eo). The asymptotic behaviour of this Lanc- 
zos coefficient as s — > 00 is also demonstrated in Figure ^ 
which can be analytically shown to be 



1 — 27r/3(s) ~ Soc exp(— 47rV3 s) , 



(14) 



where the pure number Soc ~ 1.77799. This is found 
using the leading order solution for ^(e) as ^ ^ 00 



2V6 {(^^ eo) 
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(15) 



Of great interest is the Taylor expansion of /3^(s) 
about s = 0||ri|, the "Plaquette Expansion", as this 
forms the basis of the application of the general theo- 
rems to non-intcgrable models. Using the exact defining 
equations one can perform the two reversions of series 
developments, which has the initial terms 



2.7.127 . 
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(16) 
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This coincides exactly with an early calculation ||l3j and 
an independent one using the cumulants directly, al- 
though we have calculated many more terms. To indi- 
cate the asymptotic behaviour of the coefhcients in this 
Taylor expansion, and the convergence of the series we 
have plotted the nth root of the coefhcients in Figure |^ 
versus the index. There is a clear finite limiting value in- 
dicating a finite radius of convergence for the expansion, 
and an interesting quasi-periodic cusp structure super- 
imposed. Another indication of this convergence is given 
in Figure ^ where the complex zeros of the truncated 
/3^(s) are plotted in the complex s-plane for a sequence 
of truncation orders. Clearly here the distribution of 
zeros approaches a limiting curve, with a simple shape 
about s = and at a finite distance from the origin. As 
the Lanczos Process terminates when /3„ — > at some 
nonzero n, this curve gives the breakdown of the Plaque- 
tte Expansion (the exact /3^(s) for real s > is nonzero). 
One can in fact prove this exactly from considerations of 
the analytic behaviour of e(^) as a function of complex 
^. This transformation of the ^-plane is analytic except 
on the imaginary axis > 7r/2 with a sequence of 

branch points at £,k = i{k+ ^)tt. The real part of e has a 
discontinuity across these cuts (the defined value on the 
cut is zero) , and both parts have inverse square-root sin- 
gularities when approaching the branch points from one 
direction along the imaginary axis. The closest zero of 
e(^) to the origin is relevant in determining the radius 
of convergence and this occurs a.s ^ ^ (in fact all the 
zeros lie on the imaginary axis, some at ^^) and therefore 
the value of sq is given by 
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So = - 
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(17) 



where (7 = 3^. The integrand has the following series 
development 



iye((T) = h - > 



fe=A;o + l 



1 - 



(fc+i)7r 



^(fc+i)27r2-a2j(18) 



where fco = It^ —\\ ■ Using this it is easy to show that 
2 



So 



■log n 



odd l>0 



0.2396967390143423611761660826222 



(19) 



Another relevant set of zeros are those of d0^{s)/ds 
and these are plotted in Figure I for the same sequence 
of truncations. At finite truncations a finite minima or 
a complex conjugate pair of stationary points (i.e. an 
inflection point) oi (3'^{s) always exists (this pair he close 
to the positive real axis), and these can be seen to lie 
just inside the curve where the Lanczos process breaks 
down. Finally we plot the relative error in the ground 
state energy estimates in Figure || using the truncated 



series directly in the theorem, Eq. (g), versus truncation 
order. After an initial reduction in the error to the one 
percent level, the error just continues to fluctuate about 
this level. The reason why no further improvement is 
possible with a direct application of the ground state the- 
orem is that the minima (or inflection point) is pinned 
just inside the radius of convergence, whereas it should 
tend to infinity in order to converge. 

We have found the exact solution to the Lanczos 
process for a quantum Many-body System in the ther- 
modynamic limit using results from the theory of random 
Matrix Ensembles and Orthogonal Polynomial Systems. 
Our exact method applies to any model where the cumu- 
lant generating function can be found, where it should 
be noted that in this case it was not necessary to fully 
diagonalise the problem in order to find this. Our exact 
solution provides an illustration of some general theorems 
applying to arbitrary extensive Many-body Systems and 
also demonstrate that these theorems continue to work 
as an approximate technique even at a critical point in 
the model, albeit to only a semi-quantitative level. The 
analytic Lanczos method has recently been applied to 
a related model with a quantum phase transition, the 
Ising model in a transverse field 15 1, to elucidate the 



dependence of the Lanczos Process away from, near and 
at criticality. This work was supported by the Australian 
Research Council. 
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